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Abstract. The statistical properties of the spectrum of the staggered Dirac operator in an SU(2) lattice 
gauge theory are analyzed both in the bulk of the spectrum and at the spectrum edge. Two commonly used 
statistics, the number variance and the spectral rigidity, are investigated. While the spectral fluctuations at 
the edge are suppressed to the same extent as in the bulk, the spectra are more rigid at the edge. To study 
this effect, we introduce a microscopic unfolding procedure to separate the variation of the microscopic 
spectral density from the fluctuations. For the unfolded data, the number variance shows oscillations of the 
same kind as before unfolding, and the average spectral rigidity becomes larger than the one in the bulk. 
In addition, the short-range statistics at the origin is studied. The lattice data are compared to predictions 
of chiral random-matrix theory, and agreement with the chiral Gaussian Symplectic Ensemble is found. 
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1 Introduction 



The spectrum of the Dirac operator is an important as- 
pect of nonperturbative QCD. In Euclidean space, the 
Dirac operator reads ip = i^ + g{X°-/2)^°', where g is the 
coupling constant, A° are the generators of SU(A^)-color, 
and A° are the gauge fields. Whether in the continuum 
or on the lattice, to obtain physical observables one has 
to perform an average over the ensemble of gauge field 
configurations. The eigenvalues of ip fluctuate over this 
ensemble of gauge fields. In many areas of physics, e.g., in 
classically chaotic systems, complex nuclei, and disordered 
mesoscopic systems, the study of energy level statistics 
plays a prominent role. Although the global spectral fluc- 
tuations are strongly system dependent, the local spectral 
correlations on the scale of the mean level spacing are uni- 
versal and only depend on the global symmetries of the 
system under consideration. These universal correlations 
can be obtained exactly in random-matrix theory (RMT) 
|l| . In order to describe the universal statistical properties 
of the eigenvalues of ip, one needs random-matrix models 
which take the chiral symmetry of the Dirac operator into 
account pi. In a chiral random-matrix model, the matrix 
of the Dirac operator (in a chiral basis in Euclidean space) 
is replaced by a random matrix W with appropriate sym- 
metries, 



ip 





M/t 



W 




(1) 



The average over gauge fields is replaced by an average 
over the ensemble of random matrices, and the gluonic 
part of the weight function, cxp(— Sgi), is replaced by a 
simple Gaussian distribution of the random matrix W. 
Depending on the number of colors and the representation 
of the fermions, the Dirac operator falls into one of three 
universality classes which were classified in Ref. Q and 
correspond to the three chiral Gaussian random-matrix 
ensembles. These are the chiral Gaussian orthogonal (ch- 
GOE), unitary (chGUE), and symplectic (chGSE) ensem- 
ble in which the random matrix W is either real, complex, 
or quaternion real, respectively. For a recent review of im- 
portant results, we refer to Ref. Q. 

Various quantities, in addition to the spectral corre- 
lation functions, have been employed to describe the sta- 
tistical properties of the spectrum of a given system. One 
uses the number variance and the spectral rigidity, to be 
defined below, to measure long-range correlations, whereas 
the distribution of spacings between eigenvalues is used to 
probe short-range correlations in the spectrum. A statis- 
tical analysis of the spectrum of the Dirac operator was 
first performed using data obtained by Kalkreuter in an 
SU(2) lattice gauge theory with both Wilson and Kogut- 
Susskind fermions g, and it was found that the lattice 
data were described by RMT. Recently, high-statistics lat- 
tice data obtained by Berbenni-Bitsch and Meyer were 
analyzed. Again, perfect agreement with RMT was found 
Q. In the bulk of the spectrum, the local spectral fluctu- 
ation properties follow the predictions of the conventional 
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Gaussian ensembles, which in this case are identical to 
those of the chiral ensembles ]j . 

The QCD Dirac operator has a special symmetry which 
has important implications for the edge of the spectrum. 
Since the Dirac operator anti-commutes with 75 , its eigen- 
values appear in pairs ±A leading to level repulsion at 
A = 0. Based on an analysis of Leutwyler-Smilga sum 
rules |8|, Shuryak and Verbaarschot g| conjectured that 
the so-called microscopic spectral density of the Dirac op- 
erator, defined by 



Ps{z) 



lim 



— P(—) 



(2) 



should be a universal quantity as well. Here, p{X) = 
{J2i S{X — Xi))A is the overall spectral density of the Dirac 
operator, V is the space-time volume, and S is the abso- 
lute value of the chiral condensate, respectively. According 
to the Banks-Casher formula, S = Trp{0)/V [||, the spec- 
tral density at zero virtuality is proportional to the chiral 
condensate, therefore the spectrum edge is of great inter- 
est for the understanding of the spontaneous breaking of 
chiral symmetry in QCD. The definition of Eq. (g) leads 
to a magnification of the region of small eigenvalues by a 
factor of VU, thus the microscopic spectral density also 
provides information about the approach to the thermo- 
dynamic limit. 

If ps is a universal quantity, it should also be calculable 
in chiral random-matrix theory (chRMT). This was done 
for all three chiral ensembles in the chiral limit [ [lO||ll] , [l4 . 
Several pieces of evidence have lent support to the univer- 
sality conjecture for ps. We mention here the reproduc- 
tion of Leutwyler-Smilga sum rules pC| , the microscopic 
spectral density in an instanton liquid model iQ, the 
valence quark mass dependence of the chiral condensate 
||lj,^, theoretical universality proofs with regard to the 
probability distribution of the random matrix [f6|jlj,0] , 
and finite temperature calculations |l9[ . Very recently, the 
conjecture was verified directly by lattice calculations. In 
Ref. iQ , the microscopic spectral density, the distribution 
of the smallest eigenvalue, and the microscopic two-point 
correlator were constructed from lattice data calculated 
by Berbenni-Bitsch and Meyer for an SU(2) gauge theory 
with staggered fermions in the quenched approximation 
and compared with the random-matrix predictions. The 
agreement was excellent. 

Another very interesting area of application of chiral 
random-matrix theory is the construction of schematic 
models for the chiral phase transition at finite tempera- 
ture and/or chemical potential |l|,|2|J2|||,||,|6|J23,|| . 
It should be emphasized that while chRMT yields ex- 
act results for the spectral correlations in the bulk of the 
spectrum and at the spectrum edge, the results of these 
schematic models typically are not universal. Neverthe- 
less, such models can shed light on a number of important 
issues, such as the nature of the quenched limit at finite 
chemical potential p3|. Other aspects concern the rela- 
tion with Nambu-Jona-Lasinio models pfl] and the Ua(1) 
problem [g9|. In this paper, however, we shall not discuss 
such models but concentrate on the universal properties 
of the spectrum of the Dirac operator. 



It is the goal of the present work to perform the statis- 
tical spectral analysis at the spectrum edge for the same 
lattice data as in Ref. |g^ and to compare to the statisti- 
cal properties in the bulk. Our aim is two-fold: (i) to test, 
for higher order correlations, the conjecture that the spec- 
tral statistics of the QCD Dirac operator are universal and 
described by chRMT, and (ii) to gain an intuitive under- 
standing of the statistical properties of the Dirac spectra 
at the edge (at zero virtuality). The difference between 
the usual situation (normally in the spectrum bulk) and 
the present one is the hard edge at A = which breaks the 
translational invariance of the spectrum in this region. 

In Sec. g we give, for the microscopic region, the the- 
oretical predictions for the one- and two-point spectral 
correlation functions of the chGSE and the expressions 
for the number variance and the average spectral rigidity 
applicable to the present situation. Section || is devoted to 
the analysis of lattice data with regard to the long-range 
statistics, as well as the one- and two-point correlation 
functions. We then perform, in Sec. H, a microscopic un- 
folding in order to obtain a better understanding of the 
behavior of the long-range correlations obtained in Sec. O. 
In Sec. ra we study the spacing distribution in a specific 
case. The last section is a summary and also discusses the 
effect of lattice parameters on our results. 



2 Theoretical predictions of the chGSE 

The appropriate ensemble corresponding to staggered fer- 
mions in SU(2) is the chGSE Q. A central object in 
RMT is the joint distribution function of the eigenvalues 
of the matrix in Eq. (|l|) which is obtained after diagonal- 
izing the matrix and computing the Jacobian. If W has 
N + 1/ rows and N columns, then the matrix in (mj has 
N positive eigenvalues Ai, . . . , Aat, iV negative eigenvalues 
— Al, . . . , — Atv, and i^ zero modes (we assume v > 0). The 
joint distribution function then reads § 



N 



N 



P(Ai,...,Ajv)^n^^ 



2Nf+4iy+3 -2NE^\] 



n(^'-^?)' 



k<l 



(3) 



where we have omitted a normalization factor which en- 
sures that J dXi ■ ■ ■ cIXnp = 1. Here, Nf is the number 
of massless flavors, and ly corresponds to the topological 
charge. For the quenched lattice data we consider, Nf = 
and 1^ = 0. 

We now need expressions for the scaled one- and two- 
point functions at the spectrum edge. They can be derived 
from results computed by Nagao and Forrester |12| for the 
Laguerre ensemble with the joint distribution function 



N N 

P(xi, . . . , xjv) ~ n 4"e-4^^- Y[{xk - xi)^ (4) 

j=l k<l 



by a simple transformation of variables, Ax — 2N2J^X^. 
Evidently, 4a — Nf + 2v + 1. We use the same letter P 
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to denote the joint distribution functions (0) and (0). For 
convenience, we introduce a new parameter 



^^ Nf + 2v 



(5) 



1 /.I 

2 



which will characterize the dependence of the various quan- 
tities we consider on Nf and v. For the present lattice data 
we have /i = 0. The /c-point function follows by integrating 
the joint distribution function over all but k variables, 



with 

S{x,y)^ I duu^ I dv[Jf^{2uvx)Jfj^+i{2uy) 

vJ^{2ux)J^+i{2uvy)] , (11) 
I{x,y) = / duu / dv[Jfj,{2uvx)Jfi{2uy) 



RkiXi,...,Xk) = 

Nl 



{N -ky. 



dXk+i ■ ■ ■ dAjvP(Ai , . . . , Aat) . (6) 



Jf,{2ux)Jf,{2uvy)] , (12) 



A note on notation: We denote the usual fc-point func- 
tions by Rk and the corresponding microscopic limits by 
Pk- The microscopic limit of the fc-point function is ob- 
tained by rescaling all arguments by a factor of 2NS in 
analogy with Eq. y^. Setting the lattice constant to unity, 
we can identify the volume V with the number of eigen- 
values 27V ^ V K, 2N in the large- iV limit. Including the 
proper normalization and using the result of Ref . p2[ , we 
obtain for the microscopic one-point function 



D{x,y) 



1 ,.1 

3 



duu / dvv[J^jf-i{2uvx)Jf^-f-i(2uy) 



Jfj,+i{2ux)Jfj,+i{2uvy)] . (13) 



Again, this result can be simplified further and expressed 
in terms of two single integrals |Q . We obtain 

^2(^1,^2) = fiTrzi,Trz2)d^-,d22f{7TZi,Trz2) 

+ dz^f{TTZl,TTZ2)dz2f{T^Z2,TTZl) (14) 



Ps{z) 



1 



-M[ 



/ z 



N^^oo2NS"' \2NEj 

2z^ / du'iP' / dv\J ^(2uvz)J p,j^\(2uz) 
Jo Jo 



with 



f{x,y) = ^J dtC{t,2y)--J dtC{t,2x) (15) 



-vJf,{2uz)J^,+i{2uvz)] , (7) 

where J denotes the Bessel function. Asymptotically, Ps{z) 
— » I/tt as z — > (X). To compare with the bulk properties, it 
is convenient to rescale the argument by a factor of tt so 
that it is measured exactly in terms of the local mean level 
spacing. (This factor of tt comes in through the Banks- 
Casher relation, S = Trp{0)/V.) The result in Eq. (m) can 
be simplified further and expressed in terms of a single 
integral |Q. We obtain after some algebra 

pi{z) = irpsiTTz) 



and 

^/ N ^ xJf,+i{x)Jf,{y) - yJf,{x)J^+i{y) 
x'^ — y"^ 

The derivatives of / can be expressed as 

1 f^y 

dj{x, y) = 2yC{2x, 2y) - -J^{2x) / dtJ,,{t) (17) 



and 



dxdyf{x,y) 



7T^Z[JI{27TZ) - J,,+ i(2^z)J^_i(2^z)] 



9 2 

X — y 



P2-ITZ 

-|j^(27rz) / dtJ^,{t) 



(8) 



The two-point spectral correlation function is given by 

i?2(Ai,A2) = i?i(Ai)i?i(A2) - r2(Ai,A2) , (9) 

where r2(Ai,A2) is the two-point cluster function which 
contains the non-trivial correlations. We write the mi- 
croscopic limit of Eq. (||) as ^2(^1, ^2) = Pi{zi)pi{z2) — 
T2izi,Z2). We are mainly interested in the microscopic 
limit of T2, i.e., in r2. Making use of the results of Ref. iQ, 
we obtain 



r2iz„Z2) = ^un^{^yT2 



TTZl 7rZ2 



2iVZ" 2NU 



x[{x^+y^)[2C{2x,2y) ~ J^{2x)J^{2y)] 

+xy [j^+i (2x) J^_i (22/) + J^_i {2x)J^+i (2y)] } . (18) 

We wish to study the number statistics at the spec- 
trum edge. In the study of spectral statistics, one usually 
has to unfold the empirical spectrum in order to sepa- 
rate the global variations from the local fluctuations, since 
the former are not universal and beyond the predictions 
of RMT. In the present case, the global spectral density 
near the edge is constant to a good approximation, there- 
fore no unfolding is necessary in this region. We simply 
rcscalc the energies to introduce a dimensionless variable 
S as in Eqs. (0) and dlQ). Consider a small region [0, L]. 
The average number of eigenvalues within this region is 



{2n^ziZ2)'^[S{TTZi,Trz2)S{'KZ2,nzi) 

+ I{nzi,Trz2)D{nz2,nzi)] (10) 



N{0,L) = / d\Ri{X) . 
Jo 



(19) 
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Using the scale we set above, i.e., measuring the length of 
the interval in terms of the local mean level spacing, we 
define 



and have 



S^'-^L 



N{0,S)^ / dzpi{z) 
Jo 



(20) 



(21) 



in the thermodynamic limit N ^ oo with S finite. The 
number variance is defined as S^{I) = {{N{I) — N{I))'^), 
where / is a given interval and N{I) is the number of eigen- 
values therein. In the interval [0, S], the number variance 
can be expressed as 

S^{0,S)^ dzpi{z)~ / dzidZ2T2{zi,Z2) . (22) 

Jo Jo Jo 

Note that S"^ should not be confused with S, the absolute 
value of the chiral condensate. 

Apart from the number variance, the spectral rigidity 
A3 introduced by Dyson and Mehta has played a major 
role in the study of spectral statistics. It is defined as 
the mean-square deviation of the cumulative level density 
iV(0, L) in an interval [0, L] from the best-fitting straight 
line, 



Z\3(0,L) = min- / dx\N(0,x) - Ax - B\ 

A.B L Jq 



= {N\0,x))-{NiO,x))^ 

-^{ix^L/2)N{0,xr 



(23) 



with (• • •) = -^ Jg dx ■ ■ ■ . In RMT, the averaged spec- 
tral rigidity can be expressed in terms of i?i and i?2- 
Again, we are interested in a microscopic region at the 
spectrum edge. Using the scale defined in (po|), one has 
for the chGSE in the microscopic region 



^3(0,5) 



dzpiiz) {S^z - iS^z^ + 6Sz^ - 3z'*) 



1 

2 f^ f^^ 

— / dziZi / dz2P2izi,Z2) 

^ Jo Jo 

I rS rS 
54 



-'0 



dzidz2P2izi, Z2)[S^ {Zi+ Z2) - 4:S'^ZiZ2 



+3SziZ2{zi + Z2)-3ztzi] . (24) 



For the number variance and the averaged spectral 
rigidity in an arbitrary interval [So, Sq + S] we have similar 
expressions. These theoretical predictions from the chGSE 
will be compared with lattice data in Sec. S. We will also 
compare the statistical properties at the edge with those 
in the bulk, which should be described by the GSE as 
mentioned in the introduction. Explicit expressions can 
be found in Ref. il. 



3 Spectral correlators and long-range 
statistics 

We now analyze lattice data calculated by Berbenni-Bitsch 
and Meyer. They computed complete spectra of the stag- 
gered Dirac matrix for an SU(2) gauge theory using vari- 
ous values of (3 and a number of different lattice sizes. We 
will focus on the data from an 8^ lattice with p = 2.0. 
Here, 3896 independent configurations were obtained. Be- 
cause of the spectral ergodicity property of RMT, one can 
construct the spectral correlations in the bulk with much 
fewer configurations since the ensemble average can be 
replaced by a spectral average g. However, if one is in- 
terested in the spectrum edge one has to perform an en- 
semble average, therefore a large number of independent 
configurations is needed. 

For the present spectra, we have computed the aver- 
age global spectral density and found that it can indeed 
be separated from the local fluctuations by unfolding. We 
also found that within the small interval at the spectrum 
edge we are interested in, there is no visible variation in 
the spectral density. Therefore, no unfolding is needed in 
the microscopic region. However, we will perform a "mi- 
croscopic unfolding" later, see Sec. H. In order to compare 
the data with the predictions of chRMT, only a simple 
rescaling by the local mean level spacing is done in small 
intervals starting at the spectrum edge. A typical interval 
includes approximately 20 eigenvalues. 

Figure H shows the results for the microscopic spectral 
density. We found that the agreement with the chGSE is 
very good for z < 2, quite impressive in the region 2 < 
z < 5, and getting worse as z increases. It is interesting 
to discuss what determines the domain of validity of the 
random-matrix result. The microscopic spectral density at 
the edge is essentially 



pi{z) ^ D{0)Ri{zD{0)) 



(25) 



where D{0) is the local mean level spacing at the edge. In 
order to obtain agreement between the lattice data and the 
thermodynamic limit of the microscopic spectral density 
given by Eq. (||), z has to be small so that zD{0) <^ 1. This 
is because by taking the thermodynamic limit in Eq. (uu, 
we treat the argument A of Ri (A) as an infinitesimal quan- 
tity. Thus, the random-matrix result is valid over a larger 
range if D{0) is_sman. We have D{0) = 1/VU, and the 
actual value of D{0) is determined by several factors. It 
decreases with increasing lattice volume and increasing 
chiral condensate (everything is in lattice units). The con- 
densate (at fixed lattice volume) depends on the coupling 
constant and on temperature, and it decreases with both /3 
and T. Thus, the domain of validity of the random-matrix 
result is smaller for smaller lattice volume, larger /3, and 
larger temperature. To quantify these statements would 
require a systematic study of lattice data at several values 
of V, /?, and T which is beyond the scope of this paper. 
For the present lattice, we found D{0) « 0.0059. There- 
fore, agreement with Eq. (||) is restricted to z ^ 170. In 
practice, we observe agreement only for values of z much 
lower than this upper bound. As we stressed earlier, the 
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Fig. 1. The spectral density at the edge on the scale of the 
mean level spacing calculated from the lattice data and from 
the chGSE prediction is compared (a) with that in the bulk 
and (b) with the chGUE prediction. 



volume Y has to be identified with 2A^, i.e., with twice the 
dimension of the matrices in the random-matrix model. 
Of course, one can also construct p\ for purely random 
matrices of finite dimension. In this case, the agreement 
with Eq. (|^) is quite good already for small dimension. In 
the case of the lattice data, however, the agreement with 
Eq. (^ is worse for finite V because the lattice Dirac op- 
erator has additional non-random components. We thus 
attribute the disagreement between the lattice data and 
the prediction of chRMT for z > 5 to both the validity 
of the thermodynamic limit and the non-random compo- 
nents of the lattice Dirac operator. 

From Fig. |l|, we see that the microscopic spectral den- 
sity has an oscillatory pattern with peaks distributed al- 
most periodically with period k, 1. The position of the i-th 
peak is the most probable value of the i-th eigenvalue. The 
distribution of the eigenvalues at the edge looks somewhat 
like a picket fence as does the distribution of energy lev- 
els of a harmonic oscillator. This is a consequence of the 
strong repulsion of the eigenvalues with the fixed point 



Fig. 2. The averaged staircase function at the spectrum edge 
and in the bulk. The full and the long-dashed line represent the 
chGSE and the chGUE, respectively. The short-dashed line is 
N — S, representing the averaged staircase in the spectrum 
bulk. The symbols are the data. 



z — 0. The height of a peak decays as z increases and 
asymptotically tends to 1. In Fig. |^(a), the spectral den- 
sity in the bulk (around the midpoint of the spectra) is 
also plotted for comparison. There are only slight fluctua- 
tions around the average value pbuik = 1- We also plot the 
chGUE result for the microscopic spectral density |10|, 



piW- Y^i'^oM + '^i'M] 



(26) 



for Nf = v — 0, in Fig. |^(b). It can be seen that the peaks 
of the chGUE are less pronounced and decay faster than 
those of the chGSE due to its weaker repulsion (Dyson 
index ^ = 2 vs ^9 = 4). 

In Fig. |2[ we plot the averaged number of eigenvalues in 
the interval [0, S], i.e., the averaged staircase, as a function 
of S at the edge as well as in the bulk. Again, the agree- 
ment with the chGSE, Eq. (|2l|), is good. The short-dashed 
line represents the homogeneous distribution N = S in 
the bulk. (Note that for the lattice data we consider, the 
number of eigenvalues is large enough so that for the small 
range of S shown in Fig. |2| unfolding in the bulk is not 
actually necessary.) The deviation of the chGSE from the 
straight line N = S at the edge is much larger than that 
of the chGUE as in the case of the microscopic spectral 
density. 

We now turn to the two-point correlator. In Fig. ^ we 
plot the results for P2{zi,z2) and T2{zi,Z2) as a function 
of zi for two fixed values of Z2, -^2 = 1-6 and Z2 — 2.6. One 
can see that the histogram for p2(-Zi, ^2) agrees well with 
the chGSE prediction, although the agreement is not as 
good as in Fig. H for the one-point function. For T2{zi,Z2), 
the statistics are worse than for /52 (2^1,^2), except for a 
small region around zi = Z2- This is because P2{zi, Z2) also 
includes the one-point functions which are dominant and 
have better statistics. From the definitions of p2{zi,Z2) 
and T2 (2:1,22) we can see that one needs a much larger 
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Fig. 3. The microscopic two-point correlation function 
P2{zi,Z2) and the corresponding microscopic two-point clus- 
ter function T2{zi,Z2) as functions of zi for fixed Z2. In (a) 
and (c), 22 = 1.6, in (b) and (d), Z2 — 2.6. The histograms 
represent lattice data, the dashed lines the chGSE predictions. 



number of spectra to obtain better statistics. For a given 
value of Z2 , only those configurations in which there is an 
eigenvalue in the bin around Z2 contribute to the two-point 
function. Even for a value of Z2 chosen at a peak of pi (as in 
Fig. |), no more than 1/3 of the configurations are actually 
involved in the ensemble average in the construction of 
P2(^i, Z2) and T2(zi, Z2) from the data. This is the reason 
why we only computed p2{zi, Z2) and T2{zi, Z2) for these 
specific values of Z2. For other values, one would have to 
choose a larger bin size. 

Figure shows our results for the number variance 
S'^{So, So + S) at the edge for different values of S'o com- 
pared with the corresponding chGSE predictions. We also 
calculated the number variance in the bulk for compari- 
son in Fig. H(a). No spectral average was performed in the 
bulk. From this figure, one observes the following points: 
(i) At the edge, significant systematic deviations from the 
theoretical prediction occur when S'o + S' > 4. The data 
approach the asymptotic behavior faster than the theo- 
retical curve. Again, this means that the agreement with 
the thermodynamic limit is restricted to the small re- 
gion [0,4 ~ 5], cf. the discussion after Eq. (p5|). (ii) For 
S ^ 1, the number variance at the edge increases very 
slowly, in contrast to the linear relation S"^ = S in the 
bulk. This is simply a manifestation of the suppression 
of the microscopic spectral density in this region. In fact, 
from Eq. ( |2^ ) it can be seen that S^ is dominated by the 
first term, i.e., by the average staircase function, when 
5 <C 1. The two-point correlations manifest themselves in 
U^{So,So + S) only ii So + S > 1. (ui) The overall value 
of U^{So,So + S) strongly depends on the value of So- 
The larger the value of pi{So), the larger S^{So,So + S). 
Two extreme cases are shown in Fig. |4|(a) where the entire 
curves of S'^{So,So + S) are higher resp. lower than the 
GSE curve. The values S'o = 0.63 and So = correspond 
to the first maximum and minimum of pi, respectively. 
Figures |j(b) , (c), and (d) show the results for So = 0.83, 
0.96, and 0.89 with pi{So) < 1, > 1, and = 1, respectively, 
(iv) For fixed So, the curves show strong oscillations with 
peaks appearing when pi(So + S) reaches its maxima. The 
values of S'^ for arbitrary intervals [Si, S2] at the edge are 
distributed around the GSE curve. It is known |3^ that 
the number variance refiects the "compressibility" of the 
eigenvalue "gas". Therefore, roughly speaking, the com- 
pressibility at the edge is on average the same as in the 
bulk. The case So = demands special attention. The very 
strong suppression shown in Fig. ^(a) is special because 
the interval starts at the origin which is a fixed point for 
all spectra so that fiuctuations of the eigenvalue number 
from the left hand side are prohibited. It is always harder 
to compress this one-dimensional gas on one side than on 
two sides. 

At first glance, one might attribute these features of E"^ 
to the inhomogeneity of the microscopic spectral density. 
To clarify this point, we will in the next section introduce a 
microscopic unfolding procedure to remove the oscillations 
in the microscopic spectral density, and then investigate 
the effect of the remaining two-point correlations on Z"^ 
and A3. 
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Fig. 5. The averaged spectral rigidity at the edge in intervals 
(a) [0, S\ and (b) [0.83, 0.83 + S\ as calculated from the lattice 
data is compared with chGSE (solid lines) and GSE predic- 
tions. 
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Fig. 4. The number variance in an interval [5o, 5o + S\ at the 
edge with (a) ^o = and 0.63, (b) ^o = 0.83, (c) So = 0.96, and 
(d) So = 0.89 as calculated from the lattice data is compared 
with the corresponding predictions of the chGSE (solid lines) 
and GSE (dashed lines). For comparison, the number variance 
in the bulk is computed in (a). 



We have also calculated the averaged spectral rigidity 
both at the edge for various intervals and in the bulk. In 
Fig. I we plot ^3(0, S) and ^3(0.88, 0.83 + S) compared 
with the corresponding chGSE predictions. We see that 
the region where the lattice data agree with the predic- 
tion of the chGSE is larger than in the case of the num- 
ber variance, and that the agreement is nearly perfect. 
In contrast to the case of the number variance, we find 
that the averaged spectral rigidity at the edge is always 
smaller than that in the bulk no matter what value of ^o 
is chosen, indicating that the spectrum at the edge is more 
rigid than in the bulk. Moreover, the difference between 
the edge and the bulk for the averaged spectral rigidity is 
small, compared to the significant difference in the case of 
the number variance. For S" <C 1, the behavior of ^3(0, 5') 
is dominated by the one-point function and lower than 
the straight line A-^ = S/lh in the bulk case. In addition, 
the curve shows slight convex-concave oscillation with the 
same period as the oscillations in the number variance and 
the microscopic spectral density. 
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4 Microscopic unfolding and two-point 
correlations 

As noted in the previous section, both the number vari- 
ance and the averaged spectral rigidity at the spectrum 
edge show oscillations with the same period as the micro- 
scopic spectral density. We now wish to clarify whether 
these oscillations are simply the manifestation of the in- 
homogeneity of the microscopic spectral density, or are 
intrinsic in the two-point correlations at the edge. To this 
end, we first need to separate the variation of the spectral 
density from the fluctuations on a smaller scale than the 
microscopic one. This procedure is known as unfolding. 
We therefore unfold the eigenvalues of all spectra at the 
edge using the ensemble average of the staircase function. 



unfolded chGSE 

chGSE 
GSE 




N{z, 



dzpi{z) 



(27) 



Here, Zi is the z-th eigenvalue measured in terms of the 
local mean level spacing. Since theoretical results for the 
staircase function and the microscopic spectral density are 
available, we can use them directly to unfold the data. The 
function pi{z) is given by Eq. (||). 

The microscopic unfolding procedure now works as 
follows. We introduce a variable transformation from z 
(which is on the microscopic scale, cf. Eq. ([20|)) to a new 
variable z' , 



(28) 



+ 0.4 

- o 



w 



Z ^r z' = I dtpi{t) = Oj{z) 

Jo 
We denote the inverse of u!{z) by i^{z'), 

z = n{z') . 



0.0 





unfolded ChGSE 


(b) 


ChGSE 

GSE 


o/^-°o °/A°o °°/^ 


J 



(29) 



3 
S' 



In terms of the new variable z' , the microscopic spectral 
density becomes 



dz 



p[iz')=p,iz)—^p^{Q{z')) 



1 



dz' 



:JV^=' (30) 



as it should be. The two-point correlation function be- 
comes 

1 



P2(%>^2) = P2{zi,Z2)-r7 



(31) 



dz'i/dzi ■ dz2l dz2 

P2{Q{Z[),Q{Z'2)) 

p^{n{z[))p^{Q{z'^)) ' 
and the two-point cluster function is 

pi{f2[z[))pi{fi{z'^)) 

The number variance in an interval [0, S'] on the unfolded 
scale is given by 

r'2(0,S") = S"- / / dz[dz'^T!^{z[,z'^) 
Jo Jo 

dzpi{z) - I I dzidZ2T2{zi, Z2) 

Jo Jo Jo 

= ^2(0,^2(5")). (33) 



Fig. 6. The number variance calculated from the unfolded 
lattice data at the edge in an interval (a) [0, 5*'] and in an 
interval (b) [0.75,0.75 + 5"] corresponding to that in Fig. |(b). 



But for the spectral rigidity, Z\'3(0,S") ^ Z\3(0, ^2(5")). 
Instead, 

^'3(0,^') ^ ^ 4- 1^ '^''^''^C dz'2T^{z[,z'2) 



1 



S' rS' 



"'0 



dz[dz'2r^iz[,z'2)[S'\z[+z'2) 



AS' z'iz'2 + 3S' z[z'2{z[ + Z2) — 3z[ z'2 ] 



S" 2 r^^i^') f^i 

^TP + 'FJ dziUj{zi) / dZ2T2{zi,Z2) 

15 S' Jo Jo 

+ -^ / / dzidz2T2{zi, Z2){3uj{ziyuj{z2y 

~3S'uj{zi)lli{z2)[uj{zi) +Uj{z2)] 

+AS'^uj{zi)lu{z2) - S'''[uizi)+u{z2)]} . (34) 

Figure |6| shows the number variance at the edge for the 
unfolded data in the intervals [0,5"] and [0.75,0.75 + S'] 
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Fig. 7. The number variance of the chGOE, chGUE, and 
chGSE before and after unfolding in an interval [0, S\ at the 
spectrum edge. 



corresponding to the intervals [0, 5*] and [0.83, 0.83 + S'] be- 
fore unfolding as in Fig. 3, compared with the theoretical 
prediction of Eq. ( [3^ ) and its generalization to [5*0, S'o+5'']. 
On the unfolded scale, we see that the lattice data agree 
with the theory in the region < 5" + S'q < 4. As ex- 
pected, for 5" <C 1 the number variance now returns to 
the normal case, i.e., to the straight line S''^ ^ S'. For 
S" > 1, however, the curves still show oscillations with al- 
most the same amplitude as before unfolding. Also, they 
still depend strongly on the value of S'q. The only differ- 
ence is a small shift along S'. We therefore come to the 
conclusion that this kind of oscillation is intrinsic in the 
two-point correlations, rather than a simple manifestation 
of the inhomogeneity of the microscopic spectral density. 
In Fig. 0, we compare the chGSE prediction with those of 
the chGOE and chGUE under the same conditions, i.e., 
Nf — and i^ = 0. One can also see oscillations in the ch- 
GOE and chGUE curves, although they are weaker than 
those of the chGSE. 

We also calculated the averaged spectral rigidity at 
the edge for the unfolded data and found good agreement 
with the corresponding chGSE curves, as shown in Fig. H. 
We also found that after unfolding the A3 is larger than 
that before unfolding and even larger than the GSE result 
for large S'. This might imply that the smaller rigidity 
(compared to the GSE) before unfolding is mainly a man- 
ifestation of the oscillations in the one-point function. The 
picket-fence-like distribution of pi makes the spectra very 
rigid. As before unfolding, a convex-concave oscillation is 
seen. 

In Ref. 113], Nagao and Forrester also define an un- 
folding procedure. Although their definition looks quite 
different from ours, one can show that they are essentially 
identical. At first glance, their unfolding seems to remove 
the global fluctuations of the spectral density defined in 
terms of the original scale (cf. Eq. (7.7) of Rcf. ^): 




Fig. 8. The averaged spectral rigidity calculated from the un- 
folded lattice data at the edge in an interval (a) [0, S'] and in an 
interval (b) [0.75,0.75-1-5'] corresponding to that in Fig. |(b). 



However, to have a finite X, one has to consider small x so 
that X ~ 1/A^ if A'^ ^ 00. This corresponds to our rescaling 
X ^ z = 2NUX/TT before the unfolding of Eq. (|^. We 
must emphasize that because of our microscopic unfold- 
ing, the fluctuations on the scale of the mean level spacing 
(z ~ 1) have been changed. Only on the "sub-microscopic" 
scale (z <C 1) do the fluctuations remain the same as be- 
fore unfolding. For example, the level-repulsion law on the 
scale of the mean-level spacing must be the same before 
and after unfolding. But the "long-range" statistics, like 
the number variance Z''^(5q,S^ + S"), is different from 
S'^{Sq,So -1-5), as seen in Fig. ^. Nevertheless, our results 
show that for an ensemble of spectra with a homogeneous 
spectral density, the number variance can show strong os- 
cillations. Therefore, this kind of oscillation in the number 
variance can not be attributed to the oscillations in the 
microscopic spectral density but is inherent in the two- 
point correlations. 



X = 



dx'Riix') 



(35) 
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Fig. 9. The hole probability E{0, S) (a) before and (b) af- 
ter microscopic unfolding, compared with the corresponding 
chGSE predictions. The circles are the lattice data. 



5 On the spacing distribution 

We now turn to the short-range statistics. As usual, we 
define 

E{Lo,Lo + L)= j ■■■ j d\^--- d\Np{Xi, ...,\n), (36) 



where "out" stands for [0, Lo] and [Lq + L,oo). Known 
as the "hole" probability, E{Lo,Lo + L) is the probabil- 
ity that the interval [Lq, Lq + L] is free of eigenvalues. A 
related probability density, F(Lq, Lq + L), is defined as 

F{Lo,Lo + L) = 
N I ■■■ I d\2--- d\Np{Lo + L, A2, . . . , \n) ■ (37) 



F(Lq, LQ+L)dL is the probability that the interval [Lq, Lq 
+L] is free of eigenvalues and that an eigenvalue is found 
in [Lo + L, Lq + L + dL]. The nearest-neighbor spacing 



distribution function is defined as 
P{Lq,Lo + L) = 
N{N -l)f...fdX3--- dXNp{LQ, Lo + L, A3, . . . , Xn) 

out 

(38) 

so that P{Lq,Lq + L)dLQdL is the probability that two 
eigenvalues are located in [Lo, Lq + dLo] and [Lq -I- L, Lo -I- 
L -I- dL] , respectively, and that there is no other eigenvalue 
between them. Again, we rescale L ^ S = 2NSL/tt to 
move to the spectrum edge. We then have 



L;(5o,5'o + 5) = lim E 



f ttSq ttSq ttS 



N^oo \2Ni:' 2NS 2NS 



F{Sq,Sq + S)= lim -^f(^^,^^ + ^^\ , 



(39) 

I , 
(40) 



P{Sq,Sq + S) = lim 



rP 



I ttSq ttSq ttS 



N^oo {2NEY \2NS' 2NS 2NE 

(41) 

For the special case 5o = we obtain, using the results of 
Rcf. [H, 

£;(0, S) = ^51/2 exp(-7r25V2)/_i/2(7r5) , (42) 

i^(0, S) = ^53/2 exp(-7r25V2)/3/2(7r^) , (43) 

where I denotes the modified Bessel function. (The two 
quantities are related by F(0, S) = ~dE{0, S)/dS.) How- 
ever, 

P(0, S) = (44) 

due to the repulsion between an eigenvalue and the origin 
inherent in the joint distribution function, Eq. (y). By 
definition, F(0, S) is the probability density of the smallest 
eigenvalue. 

Figure |^ shows the results for £^(0, S) calculated from 
the lattice data before and after microscopic unfolding. 
We again observe nice agreement with the corresponding 
chGSE predictions. Similar results for F(0, S) are shown 
in Fig. 1^ (Note that F(0, S) plotted in Fig. |l^(a) corre- 
sponds to the distribution of the smallest eigenvalue ob- 
tained in Ref. ^^■) From Fig. |^ one can see that, after mi- 
croscopic unfolding, the shape of E{0, S) becomes closer 
to that of the usual GSE. The eigenvalues of the Dirac 
operator always occur in pairs ±Afc. Therefore, the func- 
tion L'(0, S) may also be viewed as the nearest-neighbor 
spacing distribution at the origin if one defines 2S as the 
spacing between the pair ±Amin in units of the mean level 
spacing. Since the spacing distribution depends on corre- 
lators of all order, the nice agreement with the chGSE 
predictions seen in these two figures provides evidence 
for the universality of higher-order correlations. For the 
generic P(S'o,S'o + S) defined in Eq. @) with 5o 7^ 
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Fig. 10. The distribution of the smallest eigenvalue F{0, S) (a) 
before and (b) after microscopic unfolding, compared with the 
corresponding chGSE predictions. The histograms represent 
the lattice data. 



rigidity is due to the picket-fence-like behavior of the one- 
point function, and can be reduced by microscopic unfold- 
ing. The excellent agreement between the lattice data and 
the predictions from chiral random-matrix theory provide 
direct evidence for the universality of the one- and two- 
point spectral correlations at the edge. Our study of the 
spacing distribution demonstrates this universality also for 
higher-order correlations. 

Since our study was only done for one particular lat- 
tice size and one particular value of /3, it is in order to 
discuss the effect of these two parameters on our results. 
The random-matrix results to which we compare, in par- 
ticular Eqs. (g) and ([l4|), were derived in the thermody- 
namic limit. Thus, the agreement between lattice data and 
chRMT predictions will improve with increasing physical 
volume V, i.e., with increasing lattice size and decreas- 
ing p. We have confirmed this expectation investigating 
lattice data by Berbenni-Bitsch and Meyer obtained for 
a number of /?- values in the region 1.8 •• • 2.5 and lattice 
sizes between 6'* and W^. It should be emphasized that 
for any given value of /3 one will eventually find agree- 
ment between lattice data and chRMT, provided that the 
lattice is big enough. Of course, there are practical con- 
straints. Furthermore, chRMT cannot predict a priori how 
good the agreement will be for a given set of lattice pa- 
rameters. However, the criterion after Eq. (Bq) gives an 
estimate for the quality of the agreement. The mean level 
spacing at the spectrum edge has to be much smaller than 
1/a, where a is the lattice constant. 

After all these demonstrations of universality, one in- 
teresting question is the one of practical applications. One 
possibility is the improvement of extra pol ations to the 
thermodynamic limit, as shown in Ref. |Q. Another is- 
sue is the extrapolation to the chiral limit in the presence 
of dynamical fermions. Lattice data for such an investiga- 
tion are just becoming available, and analytical work on 
the appropriate random-matrix model is in progress. 



there is, to the best of our knowledge, no theoretical ex- 
pression available. Of course, it can be constructed from 
the lattice data. We found that for a generic value of ^o, 
the histogram (which is not shown here) is quite different 
from the usual GSE prediction. It would be interesting to 
investigate this problem further. 
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denmiiller, and T. Wilke for stimulating discussions. T.W. ac- 
knowledges the hospitality of the MPI Heidelberg. This work 
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6 Summary 

We have studied the spectral statistics of the Dirac spec- 
trum in an SU(2) gauge theory with staggered fermions, 
restricting ourselves to an 8^ lattice, with emphasis on 
the long-range statistics, i.e., the number variance and the 
spectral rigidity, at the spectrum edge near zero virtual- 
ity. Our analysis shows that while the spectra at the edge 
are more rigid than in the bulk, the fluctuations are sup- 
pressed, on average, to the same extent as in the bulk. The 
strong oscillation in S"^ is an edge effect due to the two- 
point correlations, in the sense that it cannot be removed 
by microscopic unfolding. On the other hand, the larger 
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